Low-frequency measurement of the tunneling amplitude in a flux qubit 
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We have observed signatures of resonant tunneling in an Al three-junction qubit, inductively 
coupled to a Nb LC tank circuit. The resonant properties of the tank oscillator are sensitive to the 
effective susceptibility (or inductance) of the qubit, which changes drastically as its flux states pass 
through degeneracy. The tunneling amplitude is estimated from the data. We find good agreement 
with the theoretical predictions in the regime of their validity. 

PACS numbers: 85.25.Cp, 85.25.Dq, 84.37.+q, 03.67.Lx 



Several groups, using different devices, have by now 
established that superconductors can behave as macro- 
scopic quantum objects. 1-3 These are natural candidates 
for a qubit, the building block of a quantum computer. 

Qubits are effectively two-level systems with time- 
dependent parameters. One of them is a superconduct- 
ing loop with low inductance L, including three Joseph- 
son junctions (a 3JJ qubit). 4 Its potential energy, U — 
^2j=xEjj(<frj), depends on the Josephson phase differ- 
ences <f>j across the junctions. Due to flux quantization 
Ej=i = 27T$ x /$o (with $ x the external magnetic flux 
and $o — h/2e the flux quantum), only two </>j's are in- 
dependent. 

For suitable parameters, U(<pi,</>2) has two minima 
corresponding to qubit states \E' 1 and carrying oppo- 
site supercurrents around the loop. These become degen- 
erate for $ x = i<I>o. The Coulomb energy E c (= e 2 /2C, 
with C the capacitance of junction 1) introduces quan- 
tum uncertainty in the <j>j. Hence, near degeneracy the 
system can tunnel between the two potential minima. 
(Since Ec *C Ej = Ej%, we deal with a flux qubit; 
Ec S> Ej yields a charge qubit. Coherent tunneling 
was demonstrated in both.) 

In the basis {\E ,1 ,\l/ r } and near $ x = i$ , the qubit 
can be described by the Hamiltonian 

H(t) = -e(t)a z - Aa x ; (1) 

A is the tunneling amplitude. At bias e = the two low- 
est energy levels of the qubit anticross [Fig. 1(a)], with a 
gap of 2A. Increasing e slowly enough, the qubit can adi- 
abatically transform from VP 1 to \E' r , staying in the ground 
state £L. Since dE^/d&x is the persistent loop current, 
the curvature d 2 E-/d^ is related to the qubit's suscep- 
tibility. Hence, near degeneracy the latter will have a 
peak, with a width given by |e| < A. 5 We present data 
demonstrating such behavior in an Al 3JJ qubit. 

Our technique is similar to rf-SQUID readout. 6 ' 7 The 
qubit loop is inductively coupled to a parallel resonant 
tank circuit [Fig. 1(b)]. The tank is fed a monochro- 




FIG. 1: (a) Quantum energy levels of the qubit vs external 
flux. The dashed lines represent the classical potential min- 
ima, (b) Phase qubit coupled to a tank circuit. 

matic rf signal at its resonant frequency wt- Then both 
amplitude v and phase shift x (with respect to the bias 
current 1^) of the tank voltage will strongly depend on 

(A) the shift in resonant frequency due to the change 
of the effective qubit inductance by the tank flux, and 

(B) losses caused by field-induced transitions between 
the two qubit states. Thus, the tank both applies the 
probing field to the qubit, and detects its response. 

The output signal depends on the tank's quality fac- 
tor Q. Using superconducting coil, values as high as 
Q ~ 10 3 can be obtained, leading to high readout sen- 
sitivity, e.g., in rf-SQUID magnetometers. 8 Such a tank 
can therefore be used to probe phase qubits. 9 For small L, 
the results are summarized by 5 

v = Io^tLtQ/ y/1 + (2Q0 2 , (2) 
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FIG. 2: Tank phase shift vs flux bias near degeneracy. From 
the lower to the upper curve (at / x = 0) the driving-voltage 
amplitude Vdr = IoujtLtQ takes values 0.5, 1.0, 1.5, 1.9, 2.9, 
3.5, 3.9 /iV. Inset: theoretical curves for A//i = 650 MHz, 
and I = 0.07, 0.13, 0.20, 0.26, 0.39, 0.47, 0.53 nA. 
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where f x = ^x/^o — \ , lo is the bias-current amplitude, 
and fc = M/^LLt is the tank-qubit coupling coefficient, 
with M (Lt) the mutual (tank) inductance. The ground- 
state curvature is 10 
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where X(a,g) (with g = Ej/Ec) is the conversion fac- 
tor in e = Ej\f n If lo vanishes, £ = \k 2 Ld 2 E_ /d<5>i 
becomes an external parameter accounting for the qubit 
susceptibility coupled to the tank. For finite I , this has 
to be averaged over a bias cycle < <f> < 2tt. The re- 
sulting integral (4) turns out to involve a weight cos 0, 
since the effective time-dependent coupling is (kf) 2 [cf. 
the (^-derivative of Eq. (5)], proportional to the square of 
the voltage the tank induces in the qubit. The resulting 
equations are coupled and nonlinear, but readily solved 
numerically. 

For the tank, we prepared a square-shape Nb pancake 
coil on an oxidized Si substrate. The line width of the 
20 windings was 2 /xm, with a 2 /im spacing. Predefined 
alignment marks allow placing a qubit in the center. For 
flexibility, only the coil was made lithographically; an 
external capacitance Ct is used to change wt in the range 
5-35 MHz. For the selected tank (L T « 50 nH, C T « 
470 pF), we obtained w t /2tt = 32.675 MHz and Q w 725 
from the voltage-frequency characteristic. 

The 3JJ qubit structure was manufactured out of Al 
by conventional shadow evaporation. The area of two of 
the junctions was estimated using electron microscopy 
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FIG. 3: (a) Tank phase shift vs flux bias near degeneracy 
and for Vdr = 0.5 /iV. From the lower to the upper curve 
(at f x = 0) the temperature is 10, 20, 30, 50, 75, 100, 125, 
150, 200, 250, 300, 350, 400 mK. (b) Normalized amplitude 
of tanx (circles) and tanh(A/fcBT) (line), for the A following 
from Fig. 2; the overall scale k is a fitting parameter. The 
data indicate a saturation of the effective qubit temperature 
at 30 mK. (c) Full dip width at half the maximum amplitude 
vs temperature. The horizontal line fits the low-T (< 200 mK) 
part to a constant; the sloped line represents the T 3 behavior 
observed empirically for higher T. 



as 190 x 650 nm 2 while one is smaller, so that a = 
Ejs/Eji t 2 ~ 0.8. The critical current was determined 
by measuring an rf-SQUID prepared on the same chip 7 
as I c = 2eE.j/H w 380 nA. With E c /h w 3 GHz, one 
finds g w 60 and A w —4.4. The loop area was 90 /im 2 , 
with L = 38 pH. We measured v by a three-stage cryo- 
genic amplifier, placed at « 2 K and based on commercial 
pseudomorphic high electron mobility transistors. It was 
slightly modified from the version in Ref. 12 to decrease 
its back-action on the qubit. The input- voltage noise was 
< 0.6 nV/\/Hz in the range 1-35 MHz. The noise tem- 
perature was ~ 300 mK at 32 MHz. The effective qubit 
temperature due to the amplifier's back-action should be 
considerably lower because of the small k w 2 • 10~ 2 . 

The x(/x) curves measured at various Iq and a mixing- 
chamber temperature T = 10 mK are shown in Fig. 2. 
The narrow dip at / x = directly corresponds to the 



3 



one in Eq. (6), in line with the qualitative picture below 
Eq. (1). With device parameters as above, all quanti- 
ties in Eqs. (2)~(6) are known, but A only in principle: 
its exponential sensitivity to a and especially g makes 
it notoriously hard to calculate a priori. Hence, it is 
treated as a free parameter; calculated curves for the 
best fit A/h — 650 MHz are shown in the inset. For the 
largest I the experimental and theoretical curves dis- 
agree, for the rapid change of <& x then leads to Landau- 
Zcner transitions 13 ' 14 suppressing the dip. 

The T-dependence of \ is shown in Fig. 3. For increas- 
ing T the dip's amplitude decreases while, strikingly, its 
width is unchanged [Fig. 3(c)]. Both are a simple mani- 
festation of the Hamiltonian (1) yielding (<r z ) = (e/f2) x 
tanh(ft/fc B T), 15 n = Ve 2 + A 2 . This result of equilib- 
rium statistics of course assumes that the t-dependence 
of e(t) is adiabatic. However, it does remain valid if 
the full (Liouville) evolution operator of the qubit would 
contain standard Bloch-type relaxation and dephasing 
terms (which indeed are not probed 5 ) in addition to 
the Hamiltonian dynamics (1), since the fluctuation- 
dissipation theorem guarantees that such terms do not 
affect equilibrium properties. Normalized dip amplitudes 
are shown vs T in Fig. 3(b) together with tanh(A/fcBT), 



for A/h — 650 MHz independently obtained above from 
the low-T width. The good agreement strongly sup- 
ports our interpretation, and is consistent with A be- 
ing T-independent in the relevant range. 16 Of course, for 
higher T the dip will wash out; we observe a width cx T 3 
above a crossover temperature « 225 mK. For T of this 
order, deviations from the two-state model can be ex- 
pected, especially for / x ^ 0. This behavior outside the 
qubit regime has not been pursued. 

In conclusion, we have observed resonant tunneling in 
a macroscopic superconducting system, containing an Al 
flux qubit and a Nb tank circuit. The latter played dual 
control and readout roles. The impedance readout tech- 
nique allows direct characterization of some of the qubit 's 
quantum properties, without using spectroscopy. 2,3 In a 
range 50^200 mK, the effective qubit temperature has 
been verified [Fig. 3(b)] to be the same as the mixing 
chamber's (after A has been determined at low T), which 
is often difficult to confirm independently. 
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Using this simple tanh law also for the net qubit response 
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tion, flux derivatives such as in Eq. (6) become adiabatic 
rather than isothermal. For finite T and v, the problem is 
essentially nonequilibirum and outside our scope. 
S.-X. Li, Y. Yu, Y. Zhang, W. Qiu, S. Han, and Z. Wang, 
Phys. Rev. Lett. 89, 098301 (2002) and refenerces therein. 



